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Five Point Mass Gravitational Lenses in a Rhombus as a Soluble 
Model Giving the Maximum Number of Images 

Hideki Asada 

Faculty of Science and Technology, Hirosaki University, Hirosaki 036-8561, Japan 

As an extension of four point mass lenses at the vertices of a rhombus, we present five 
point mass lenses at the center and vertices of a diamond, which constitute, for a source 
behind the center, a soluble model giving expressions of all the image positions (with the 
maximum number of images as twenty for five lenses). For a source near the center, all the 
image positions are obtained in the linear approximation. 

§1. Introduction 

Since the discovery of the lensed quasar Q0957+561 in 1979 (Ref. 1), the grav- 
itational lensing has been used to find distant objects and determine cosmological 
parameters such as the cosmological constant 2 ) and inhomogeneities (e.g., Ref. 3) 
and astrophysical ones such as the mass of MACHOs (e.g., Ref. 4). 

In spite of its importance, there are unsolved problems in the gravitational lens 
theory because of its nonlinear nature. One is that roots for the so-called lens 
equation (which correspond to image positions) are unknown in analytic forms, and 
numerical roots can be obtained, except for very few cases. Very recently, pertur- 
bative methods have been provided for obtaining the roots that can be expressed as 
Maclaurin series in mass ratios. 5 ) Positions of extra images that cannot be expressed 
in Maclaurin series are still unknown as functions of lens parameters. 

Even two point masses are too complicated to treat by hand, since the lens 
equation is in fifth order in real variables 6 ) as well as complex ones. 7 ) Fifth order in 
complex variables could lead to higher orders in real ones but this is not the case. 
Only for the special position of the source, the binary lens equation becomes exactly 
solvable. 8 )' 9 ) Until recently, it has been unknown how many the maximum number of 
images for N point mass lenses is. Eventually, 10 )~ 12 ) Rhie showed that it is 5(iV — 1). 
The maxima are achieved indeed by a binary as five. One important example is four 
lenses at the vertices of a rhombus but with numerical methods. 10 ) This example 
was considered as an extension of the earlier work by Mao, Petters and Witt (Ref. 
13), where they investigated properties of point mass lenses on a regular polygon 
with N > 3, though for this model the number of images is not more than 3N + 1 
(less than 5 (A — 1)) and thus does not provide the maximum number. Rhie (2003) 
first considered a model of N equal masses on a regular polygon and next added a 
small mass at the center in order to prove that the maximum number of images is 
no less than 5(N — 1) without giving expressions for the image positions. 11 ) 

The purpose of this brief paper is to present certain soluble models for gravita- 
tional lenses. First, we shall reexamine, especially in exact algebraic manners, four 
point mass gravitational lenses at the vertices of a rhombus. Next, as an extension, 
we investigate five point mass lenses at the center and vertices of a diamond, which 
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constitute, for a source located behind the center, a soluble model giving expres- 
sions of all the image positions (with the maximum number of images as twenty for 
N = 5). In this sense, the present paper can be considered as an extension of the 
previous papers. 10 - 1 ' 11 ) 



§2. Soluble models giving the maximum number of images 

The lens equation, expressed in complex quantities, 14 )' 15 ^ is given for N point 
mass lenses by 7 ) 

N 

p =i z e p 

where the source, image and lens positions are denoted as w = f3 x + if3 y , z = 6 x + i8 y , 
£p = e px + ie py , respectively, the mass ratio is denoted as m p , and the asterisk * 
means the complex conjugate. Here, the angle is normalized by the Einstein ring 
radius defined by the total mass. 

2.1. Four lenses at the vertices of a rhombus 

Let equal masses (rrii = 1/4) located at the vertices of a rhombus (See the next 
subsection for inequal mass cases), where its center is chosen as the center of the 
coordinates and the location of each mass (in the anti-clockwise direction) is denoted 
as e\ = k, €2 = it, £3 = —k, €4 = — it for k,t G R. Without loss of generality, we 
assume k > t. 

Let us put the source at the center (i.e. the maximally symmetric location). 
Then, the equation is rewritten as 

z(z* 2 - k 2 )(z* 2 + 1 2 ) = z* (V 2 + h^t 2 - k 2 )"j . (2-2) 

We employ the polar coordinates as z = rexp(icp). 

An immediate solution for Eq. (|2-2p is r = 0. We shall seek other solutions with 
r 7^ 0. Eq. f|2-2|) is reexpressed as 

r 4 - r 2 = + B 4 e 2 ^, (2-3) 

where real quantities are defined as 

A 4 = k 2 t 2 , (2-4) 

B, = (t 2 - k 2 ) Q - r 2 ) . (2-5) 

Here, the L. H. S. of Eq. (|2-3p is real and hence the R. H. S. must be real. We thus 
have three cases as exp(2i<f>) = ±1 or j44[exp(2i</>) + exp(— 2icj>)] + B4 = 0. 
Case 1 . exp(2i<p) = 1 ((f) = or it): 
The lens equation becomes 

r 4 + { f _ k 2_ 1)r 2 + h k 2 _ e) _ f,2 ( 2 = Q; (2 . 6) 
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where the coefficient of r 2 is always negative because of k > I. This equation 
quadratic in r 2 can be solved immediately. According to Descartes' rule of signs 
(e.g., Ref. 16), which states that the number of positive roots either equals to that 
of sign changes in coefficients of a polynomial or less than it by a multiple of two, 
we have the maximum number of positive roots as 2 if {k 2 — £ 2 )/2 > k 2 £ 2 , whereas 
we have only one positive root otherwise. It is worthwhile to mention that the 
number of the positive roots can be precisely counted by using an explicit form of 
the roots, though it needs straightforward but lengthy calculations (especially for a 
cubic equation discussed later). 
Case 2 . exp(2#) = -1 (0 = vr/2 or 3vr/2): 
The lens equation becomes 



where the sum of the constant terms is always negative for k > £. This equation 
for r 2 can be solved immediately. Regardless of the sign of the coefficient of r 2 , the 
number of positive roots is one according to Descartes' rule of signs. 
Case 3 . A 4 [exp(2i(j)) + exp(-2ic/>)] + B A = 0: 

Given r, this provides four roots of <j) (mod 2ir). For this case, Eq. (|2 3|) is reduced 
to 



If and only if k£ < 1/2, we have two positive roots. 

Therefore, we have the maximum number of images as 1+2 x 2+2 x 1+4 x 2 = 15 
if the above conditions are satisfied. The number is exactly the maximum number 
of images for iV = 4. 

A rhombus with right angles is a square, for which there is a change in the Case 
1. The sum of the constant terms becomes negative. According to Descartes' rule 
of signs, the number of positive roots is necessarily one. As a result, the maximum 
number of images does not exceed thirteen for a square model. 

Figure [T] shows fifteen image positions due to four point masses at the vertices 
of a rhombus. 

2.2. Five lenses at the center and vertices of a rhombus 

Let five masses located at the center and vertices of a rhombus, where the 
location of each mass is denoted as e\ = k, €2 = i£, £3 = —k, £4 = —i£, £5 = 
for k,£ G R. Without loss of generality, we assume k > £ again. By respecting 
symmetries regarding the principal axes of the rhombus, we take mass ratios as 
mi = 7713 = /u/2, ?7i2 = nii = u/2 and 771,5 = 1 — I 1 ~ f for < < 1, < f < 1 and 
fi + v < 1. 

Let the source located at the center such that the equation can be rewritten as 



zz*(z* 2 - k 2 ){z* 2 + £ 2 ) = z* 4 + [(1 - u)£ 2 - (1 - p)k 2 )]z* 2 - (1 - n - v)k 2 £ 2 . (2-9) 




(2-7) 



r 2 + k 2 £ 2 = 0. 



(2-8) 



We employ the polar coordinates as z = rexp(i(p). Importantly, r = is no more 
root for Eq. (|2-9p . since this position is occupied by a lens mass. Henceforth, we 
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Fig. 1. Four point mass lenses at the vertices of a rhombus and fifteen images: The source is 
located at the center. The filled squares and disks denote the positions of the lenses and images, 
respectively. We assume equal masses, and k = 0.9 and £ = 0.5. 



shall seek roots for r / 0. Eq. (|2-9|) is reexpressed as 

r 6 - r 4 = A 5 e^ + B 5 e 2i ^, (2-10) 

where real quantities are defined as 

A 5 = k 2 f[r 2 -(l-fi-u)}, (2-11) 
B 5 = [(1 - v)£ 2 - (1 - fi)k 2 ] r 2 - (£ 2 - k 2 )r 4 . (2-12) 

This has a structure very similar to Eq. (|2-3p . The L. H. S. of Eq. (|2-10p is real 
and hence the R. H. S. must be real. We thus have three cases as exp(2i(f>) = ±1 or 
j45[exp(2i^) + exp(— 2i(p)] + B$ = 0, which are investigated separately below. 
Case 1. exp(2ic/>) = 1 (eft = or tt): 
The lens equation becomes 

r 6 + {£ 2 - k 2 - l)r 4 - [(1 - v)i 2 - (1 - fi)k 2 + k 2 £ 2 ] r 2 + (l-fi- u)k 2 l 2 = 0, (2-13) 

where the coefficient of r 4 is always negative because of k > £ and the last constant 
term is positive for fj, + u < 1. Roots for this equation cubic in r 2 are given by Car- 
dano's method or alternatives. 16 ^ By Descartes' rule of signs, we have the maximum 
number of positive roots as 2, regardless of the sign of the coefficient of r 2 . 
Case 2. exp(2^) = -1 = vr/2 or 3vr/2): 
The lens equation becomes 

r 6 - (£ 2 -k 2 + l)r 4 + [(1 - v)£ 2 - (1 - fi)k 2 - k 2 £ 2 } r 2 + (l-fi- v)k 2 £ 2 = 0, (2-14) 
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Fig. 2. Five point mass lenses (filled squares) at the center and vertices of a rhombus and twenty 
images (filled disks): The source is located at the center. Here, mi = = m?, — = 6/25, 
m 5 = 1/25, and k = 0.7 and I = 0.6. 



where the last constant term is positive for \i + v < 1. This equation cubic in r 2 
is solved by Cardano's method. 16 ^ By Descartes' rule of signs, we have no positive 
root if both the coefficients of r 4 and r 2 are positive. Otherwise, it is possible that 
the equation admits two positive roots. 
Case 3 . A 5 [exp(2i(j)) + exp(-2z0)] + B 5 = 0: 

Given r, this case provides four roots of 4> (mod 2ir). Eq. (12-lOh is reduced to 

r 6 - r 4 + k 2 £ 2 r 2 - (1 - fi - v)k 2 i 2 = 0. (2-15) 

By Descartes' rule of signs, we have either three or one positive root. 

As a result, the number of images can reach 2x2 + 2x2 + 4x3 = 20 if 
the conditions mentioned above are satisfied. The number twenty is exactly the 
maximum number of images for N = 5. 

Figure [2] shows twenty image positions due to five point masses at the center 
and vertices of a rhombus. 

2.3. Approximate roots for small \w\ + 1 

Here, we comment cases. For arbitrary w, the above models are not to 

be solvable by hand. For small w, however, one can construct approximate roots as 
follows. For the general form of the lens equation as w = f(z, z*), they are expressed 
in a form as 

z = z + aw + bw* + 0(H 2 ), (2-16) 
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where zq denotes a root for w 



For our specific case as 



we obtain 



(discussed above), and a and b are denned as 
1 



1 



dz* 



>0,2o) 



2 ■ 



f(z,z*) 
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rrir. 



(2-17) 



(2-18) 



(2-19) 
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(2-20) 



(2-21) 



It should be noted that these expressions are not always valid in the entire 
complex plane for w, because we assume a smooth one-to-one map between w = 
and w 7^ cases. This correspondence is singular when 1/a = 0. This is consistent 
with a fact that Eq. (|2-20p is a magnification factor for an image at zq (i.e. the 
inverse of the Jacobian of the lens mapping). Figures 1 and 2 show non- merging 
images. Hence, the factor a seems regular at w = 0. Therefore, properties of the 
images for small \w\ may be similar to those for w = 0. Further investigations of 
large \w\ are necessary for deeply understanding of the present lens model, especially 
image properties such as image shapes near critical curves. Finally, let us mention 
a single point-mass limit as k,£ — » 0, for which Eqs. f|2- 3j) and (|2- 10|) reproduce the 
Einstein ring (r = 1). 



§3. Conclusion 



As an extension of four point mass lenses at the vertices of a rhombus, we have 
presented five point mass lenses at the center and vertices of a diamond, which 
constitute, for a source located behind the center, a soluble model giving expressions 
of all the image positions (with the maximum number of images as twenty for five 
lenses). For a source near the center, all the image positions have been obtained in 
the linear approximation. It is straightforward to extend this result to higher order 
approximations . 

The present model, though it is not realistic in astronomy, seems useful for 
theoretical investigations of lensing mostly because all the image positions (and thus 
magnification factors) are given in terms of elementary functions. Comparing with 
a binary model N = 2, a rhombus model with N = 4 or 5 has an advantage in the 
sense that one can control internal shear effects by choosing a ratio between two 
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parameters k and £ in addition to mass ratios. Further investigations are a topic of 
future work. 

Regular polygons cannot realize 5(N — 1) images. On the other hand, four lenses 
on a rhombus and five ones in a rhombus can realize them. It is thus suggested that 
shear effects are important for the increase in image numbers. A squashed polygon 
may be an interesting configuration for lensing investigations. 
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